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ABSTRACT 



Pl,: 

6 

Gravitational lensing of high-redshift supernovae is potentially an important 
source of uncertainty when deriving cosmological parameters from the measured 



brightness of Type la supernovae, especially in deep surveys with scarce statis- 
tics. Photometric and spectroscopic measurements of foreground galaxies along 
the lines-of-sight of 33 supernovae discovered with the Hubble Space Telescope, 
both core-collapse and Type la, are used to model the magnification probability 
distributions of the sources. Modelling galaxy halos with SIS or NFW-profiles 
and using M/L scaling laws provided by the Faber- Jackson and Tully-Fisher 
relations, we find clear evidence for supernovae with lensing (de) magnification. 
However, the magnification distribution of the Type la supernovae used to deter- 
mine cosmological distances matches very well the expectations for an unbiased 
sample, i.e. their mean magnification factor is consistent with unity. Our results 
show that the lensing distortions of the supernova brightness can be well under- 
stood for the GOODS sample and that correcting for this effect has a negligible 
impact on the derived cosmological parameters. 



Subject headings: gravitational lensing - supernovae: general 
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1. 



INTRODUCTION 



Having established the existence of dark energy using Type la supernovae (SNIa) (Riess 
et al. 1998; Perlmutter et al. 1999; Knop et al. 2003; Tonry et al. 2003), ongoing and planned 
supernova (SN) surveys are reaching sufficient sensitivity to explore the nature of the energy 
component of the Universe driving its accelerated expansion. One such project is the Great 
Observatories Origins Deep Survey (GOODS Giavalisco et al. 2004) supernova survey (Riess 
et al. 2004; Strolger et al. 2004), aimed at breaking the degeneracy in the cosmological 
parameters by expanding the redshift range of the studied SNe, as suggested by Goobar & 
Perlmutter (1995). However, measuring distances to z>l SNe poses additional difficulties. 
The brightness dispersion caused by gravitational lensing in the inhomogeneous Universe is 
comparable to the intrinsic spread in SNIa luminosities. Magnification of SNe due to lensing 
is a systematic uncertainty that to some extent can be cured by statistics since the mean 
magnification of a large number of sources is expected to be unity relative to an homogeneous 
universe. However, at high redshifts where statistics are currently poor and lensing effects 
potentially large, magnification bias, i.e. the preferential detection of magnified sources, could 
affect the estimates of cosmological parameters as well as the measurements of SN rates. In 
order to avoid any potential bias and reduce the scatter in e.g. the Hubble diagram, it is 
desirable to correct the measured brightness of individual sources for their gravitational 
lensing. In this paper, we use a technique described in detail in an accompanying paper 
(Gunnarsson et al. 2005) to compute the magnification for a sample of SNe observed in the 



GOODS-fields, the Hubble Deep Field North (HDFN) and the Chandra Deep Field South 



The outline of the paper is as follows. §2 contains a brief description of our method. 
In §3, the GOODS-fields and the SNe are presented, while we in §4 discuss systematic and 
statistic errors. Finally, our results are presented in §5 and discussed in §6. Throughout 
the paper, we use natural units, where c = G = 1. We use a Hubble parameter of H = 
70 kms^Mpc" 1 , a matter density of f2 M = 0.3, and a dark energy density of VI \ = 0.7. If no 
explicit redshift dependence is shown, all quantities are given with present values. Quoted 
magnitudes are Vega normalized. 



(CDFS). 



2. 



METHOD 



2.1. 



Distribution of Matter in the Universe 



The deflection and magnification of light from distant sources depend on the matter 
distribution in the Universe. Gravitational lensing effects are caused by inhomogeneities in 
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the Universe and we assume that these effects are dominated by dark matter galaxy halos. 

To account for all matter in a consistent way throughout the redshift range considered 
in this paper (z < 2), all "unobserved" matter is put into a smoothly distributed compo- 
nent. In other words, we enforce a global self-consistency on our cosmological model, where 
the fraction of the matter density not associated with observed galaxies in the data-set is 
characterized by the smoothness parameter r](z). If rj = 1, all matter is smoothly distributed 
and conversely if rj = 0, all matter is located in galaxy halos (clumps). 

Since clumps and smoothly distributed matter focus light differently, the value of the 
7] parameter affects the angular diameter distances used in computing the gravitational lens 
effects. This is discussed in more detail in e.g. Kayser et al. (1997). 

If the evolution of the total matter density in the Universe is assumed to be known, 
Pm(z) = Pm(0)(l + z) 3 , the smoothness parameter at different redshifts can be computed as 

V(Z) = l-^ y (1) 

where p g (z) is the matter density in dark matter halos. Observed galaxy luminosities can be 
translated into halo masses using the relations presented in §2.2. The density of dark matter 
residing in halos can then be estimated in redshift bins by dividing the total mass of halos 
with the corresponding volume of the bin. 

In this investigation, we use observational data obtained in the two GOODS fields, 
CDFS and HDFN (Giavalisco et al. (2004); Capak et al. (2004)). The magnitude limit 
adopted is / ~ 24.5. The redshift dependence of the r] parameter for HDFN and CDFS is 
plotted in Figure 1. The smoothness parameter increases with redshift since with a fixed 
apparent magnitude limit, the depth in absolute magnitude becomes shallower at higher 
redshift and consequently, a decreasing fraction of the matter is accounted for by observed 
objects. 



2.2. Galaxy Halos 

We assume galaxy halos to be spherically symmetric and consider two different halo 
models: Singular Isothermal Spheres (SIS) and the profile of Navarro, Frenk, and White 
(NFW Navarro et al. 1997). The density profile of a SIS 

2 

flas(r) = 2^2, (2) 
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is characterized by the line-of-sight velocity dispersion a of the galaxy. The NFW density 
profile, 

pnfw(j) = , , w.. s , — r^2> ( 3 ) 

{r/r s )(l + r/r s y 

is dependent upon the scale radius r s where approximately Pnfw ^ r ~ 2 an d p s the density 
at r ~ r s /2. For the total mass of a halo described by any of these two profiles to be finite, 
the halo must be truncated at some radius r\. We have chosen to truncate the halos at r2oo, 
defined as the radius inside which the mean mass density is 200 times the present critical 
density. The total mass of a halo is thus rri2oo, the mass enclosed within r 2 oo- For a SIS halo, 
m 2 oo can be obtained from the velocity dispersion 

4 = % (4) 

Properties of NFW halos are completely specified by TO200, since r s and p s can be obtained 
numerically from m2oo (Navarro et al. 1997). Furthermore, we assume that = mfoo- 

We estimate the velocity dispersion of each galaxy using absolute magnitudes Mb de- 
rived from observations, combined with empirical Faber- Jackson (F-J) and Tully-Fisher (T- 
F) relations for ellipticals and spirals, respectively. Ellipticals are defined as objects whose 
observed Spectral Energy Distributions (SED) are best fitted by an early type spectral tem- 
plate. Spirals are objects best-fitted by spiral or later type template SEDs. For ellipticals, 
we use the following expression for velocity dispersion derived in Mitchell et al. (2005) 

log 10 a = -0.091(M B - 4.74 + 0.85/), (5) 

where we use z' — z for redshifts z < 1 and z' — 1 for z > 1. This redshift dependence 
accounts for the general brightening of the stellar population with look-back time. At high 
redshift (z > 1), where this evolution is not well known, we assume a non-evolving brightness. 
We let the error in the derived relation be represented by the observed scatter in the Sloan 
Digital Sky Survey (SDSS) measurements (Sheth et al. 2003) 

rms(log 10 a) = 0.079[1 + 0.17(M B + 19.705 + 0.85/)]. (6) 

We use the T-F relation derived by Pierce & TuHy (1992), with correction for red- 
shift dependence calculated by Bohm et al. (2004), to derive the rotation velocity for the 
spiral/late- type population, 

log 10 Knax = -0.134(M B + 3.61 + 1.22*0, (7) 

where V max is the maximum rotation velocity for the galaxy. The observed scatter in the 
absolute magnitude around this relation is rms(M B ) = 0.41 (Pierce & Tully 1992), corre- 
sponding to 

rms(log 10 Knax) = 0.06. (8) 
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We finally convert the rotation velocity in the spiral galaxies to a velocity dispersion using 



We follow the method developed in the accompanying paper by Gunnarsson et al. (2005) 
to compute magnification of SNe using a substantially modified version of a publicly avail- 
able Fortran code, Q-LET, which utilizes the so-called multiple lens plane algorithm. The 
method is only briefly outlined here and we refer the reader to Gunnarsson et al. (2005) and 
Gunnarsson (2004) for more details. 

The lens equation for multiple lens planes relates the observed and intrinsic positions of 
a source by tracing a light-ray through planes situated at each lens redshift upon which each 
lens' mass distribution is projected. From the image plane and in every consecutive plane, 
the deflection angle is computed and the ray is recursively followed to the source plane. 
Using the Jacobian determinant of the lens equation, the magnification can also be found. 
Whereas the deflection angle only depends on the mass within the rays impact radius on 
the lens (for circularly symmetric mass density projections), the magnification also depends 
on the surface mass density at this radius. Both quantities are therefore dependent on the 
mass and density profile of the lens halo. Furthermore, the distances (which depend on the 
rj parameter) to the planes involved are important in the calculations, affecting both the 
deflection angle and magnification. 

The magnification factor ji' obtained for a specific model universe is given relative to 
a universe with the same amount of smoothly distributed matter [i.e. the same smoothness 
parameter r](z)} but with all matter in clumps infinitely far from the line-of-sight (Dyer & 
Roeder 1973) implying //' > 1 (for primary images). 

In the following, we will we present our results in terms of the magnification jj, relative 
to a universe with all matter distributed homogeneously (i] = 1). The magnifications are 
related by 



where D<r and Ds are angular diameter distances to the source calculated using 77 = 1 
(or the filled-beam approximation) and f](z), respectively (Gunnarsson et al. 2005). Due to 
flux conservation, the mean value of \x for random source positions is unity (Schneider et al. 



0- = V max /V2. 



2.3. Gravitational Lensing 




(9) 



1992). 
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3. DATA SETS 

We have examined a sample of high-redshift supernovae detected within GOODS. Of 
the total 42 SNe detected during the survey (Riess et al. 2004; Strolger et al. 2004), we 
analyze 32 SNe, divided into 19 Type la SNe and 13 core-collapse SNe. The ten SNe that 
we do not include either lack a determination of, or have uncertain redshifts, or are outside 
the area in which we have high precision photometric redshifts for the foreground galaxies. 
We also include the previously detected high-redshift SN 1997ff in our sample (Riess et al. 
2001). 

Photometric redshifts for all foreground galaxies are calculated using photometry from 
the GOODS CDFS and HDFN data sets. For CDFS we include HST ACS BViz and VLT 
ISAAC JHK S photometry (Giavalisco et al. 2004), while for the HDFN, we use KPNO- 
4m MOSAIC [/-band and SUBARU 8.2m SuprimeCam BVRIZ-h&nd data (Capak et al. 
2004), together with KPNO-4m FLAMINGOS JfT s -band data obtained in March 2003. We 
use a version of the template fitting method to derive photometric redshifts as described in 
Dahlen et al. (2005). The photometric redshift code calculates for each object the best-fitting 
redshift, the redshift probability distribution and the best-fitting spectral type. 

The accuracy of the photometric redshift code is tested by comparing with available 
spectroscopic redshifts. For GOODS HDFN, we use 848 spectroscopic redshifts taken from 
the Team Keck Treasure Redshift Survey 1 and find A z = (|z p hot — %>cc|/(l + %>cc)) ~ 0.08 af- 
ter excluding 2.4% outliers with A z > 0.3. For GOODS CDFS, we use 568 spectroscopic 
redshifts from the ESO/GOODS-CDFS spectroscopy master catalog 2 and find A z ~ 0.08 af- 
ter excluding 3.7% outliers. Both the CDFS and HDFN photometric redshift catalogs are 
complete to / ~ 24.5, i.e. we are able to determine photometric redshifts for all objects to 
this limit. At fainter magnitudes, objects start to drop out from an increasing number of fil- 
ters, making photometric redshift determination more uncertain than the quoted accuracy. 
Therefore, we only include objects to this limit. For objects with available spectroscopic 
redshifts, we replace the photometric redshifts with these. 

Rest-frame absolute magnitudes and colors are derived using the recipe in Dahlen et al. 
(2005). In summary, the absolute magnitude in, e.g., the 5-band, is calculated using the two 
observed bands that encompass the rest-frame 5-band at the given redshift. Each observed 
band is K-corrected to the effective wavelength of the 5-band using the spectral shape of 
the best-fitting template SED. The final magnitude is thereafter calculated by interpolating 



■"■http : //www2 . keck .hawaii . edu/realpublic/science/tksurvey/data43roducts/data_products . php 
2 http : //www. eso . org/ science/goods/spectroscopy/CDFS_Mastercat/ 
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between the two magnitudes, giving more weight to the filter that observe closest to the 
rest-frame £>-band, and subtracting the distance modulus. 

4. ERROR ESTIMATION 

The accuracy of our results depend on the validity of the assumption that unobserved 
matter can be treated as smoothly distributed and that our modeling of galaxy halos is 
correct. According to simulations in Gunnarsson et al. (2005) a magnitude limit of / = 25, 
approximately corresponding to the magnitude limit of the GOODS, introduce negligible 
errors in \i compared to these uncertainties. 

The largest uncertainties involved in the modeling of galaxy halos, apart from the choice 
of halo model, emerge from uncertainties in galaxy redshifts and scatter in the Faber- Jackson 
and Tully-Fisher relations. Monte-Carlo simulations were used to obtain the errors in the 
estimated magnifications due to these uncertainties. New galaxy catalogs were simulated 
based on the GOODS-fields, where the position of each galaxy was kept fixed, but redshift 
and hence absolute magnitudes were varied. In the case of a photometrically measured red- 
shift, the simulated redshift was drawn from the probability distribution of the photometric 
redshift. For galaxies with spectroscopic redshifts, the simulated redshift was drawn from a 
normal distribution with standard deviation a z = 0.01. Depending on the redshift probabil- 
ity distribution of the galaxies, foreground galaxies of a SN sometimes become background 
galaxies in the simulated catalogs and vice versa. The effects of galaxies popping in and out 
of SN lines-of-sight can be very dramatic and therefore introduce large uncertainties. When 
the velocity dispersion of a simulated galaxy was computed, the scatter in the Faber- Jackson 
and Tully-Fisher relations were taken into account. The m2oo parameter of a galaxy is, ac- 
cording to equation (4), proportional to a 3 and therefore, the scatter in the Faber- Jackson 
and Tully-Fisher relations, expressed by equation (6) and (8), were translated into a scatter 
in a 3 , or equivalently, a scatter in m 2 oo- Simulated velocity dispersions were drawn from a 
normal distribution of <r 3 -values with the translated scatter as standard deviation. For each 
simulated catalog, i](z) was computed and any negative parts of i](z) were set to zero. The 
magnification factor /x of each SN was computed for 500 simulated galaxy catalogs and the 
resulting Probability Distribution Functions (PDF) were used to estimate the uncertainties 
in the SN magnifications. 
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5. RESULTS FOR 33 SNE IN THE GOODS-FIELDS 

Figure 2 shows the foreground galaxies within 20" from the four SNIa most affected by 
gravitational lensing in the GOODS-fields: SNe 2002fx, 2003az, 1997ff, and 2003es. Redshifts 
and masses of the galaxies are indicated in the figure. In Figure 3 the magnification PDFs 
of these SNe, computed for two different halo models, are presented. The upper and lower 
panels show examples of de-magnified and magnified SNe, respectively. In general, highly 
magnified SNe have broader PDFs than moderately magnified or de-magnified ones, since 
they are more model sensitive, and thus the estimated errors increase with magnification. 
As can be seen from Figure 2, the magnified SNe 1997ff and 2003es have more crowded 
lines-of-sight than the de-magnified SNe 2002fx and 2003az. 

Collected results for the 20 GOODS SNIa are presented in Table 1. The table contains 
the magnification of the SNe and the 68% and 95% confidence intervals, computed for both 
SIS and NFW halo models, obtained by Monte-Carlo simulations. Included in the table are 
also the number of foreground galaxies within 60" to each SN. In Figure 4 the results are 
presented graphically, showing 95% confidence levels. Table 2 and Figure 5 show the same 
results for the 13 core-collapse SNe. We see that the difference between PDFs computed 
assuming SIS and NFW halos are small for most SNe but noticeable for some, like SN 2003es 
(see Figure 3). However, the PDFs computed for different halo models always overlap and 
we conclude that our method is fairly insensitive to this uncertainty. An important fact is 
that the SNe which we have found to be significantly magnified or de-magnified have errors 
smaller than the estimated magnification factor and it should thus be possible to correct for 
this gravitational lensing in a meaningful way. 



5.1. Correlation between Redshift and Magnification 

In Figure 6 the magnification, computed assuming NFW halos, of all 33 SNe is plotted 
vs. redshift. Type la and core-collapse SNe are indicated by filled circles and squares, 
respectively. The magnification scatter obviously increases with redshift for the two samples 
individually, as well as for the joint sample. Since the difference between Df 1 and 
increases with redshift, we also expect the magnification /x to be anticorrelated with redshift 
for lines-of-sight with few lensing galaxies where // ~ 1, see equation (9). This trend is 
readily seen in Figure 6. 

The scatter of magnifications expected from simulations (Gunnarsson et al. 2005) is also 
indicated in the figure. The most likely magnification as a function of redshift is indicated 
by the dashed line. Light and dark gray shaded areas represent the 68% and 95% confidence 
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level of the simulated distributions, respectively. The magnification of GOODS supernovae 
follow the redshift trend expected from simulations and the scatter is in agreement with the 
simulated data. We can thus not find any evidence for any unexpected magnification bias 
with redshift. 



5.2. Distributions of Magnifications 

The distribution of magnifications (for the NFW case) of our GOODS SNe is shown 
in Figure 7, a histogram based on the joint sample of both Type la and core-collapse SNe. 
The distribution peaks at a value slightly lower than unity and has a tail toward large 
magnifications. Moreover, the mean of the distribution is close to unity. Since the field size 
is finite, implying that not all lenses are included, and the fact that any multiple images 
are left out, we do expect the mean value of the magnification to be slightly lower than 
unity (see also Gunnarsson et al. 2005). The figure also shows, in shaded gray, the expected 
distribution of magnifications for our sample of SNe obtained by simulations. Although the 
SN sample is small, the agreement between GOODS SNe and simulated data is excellent. 

We have also compared the distribution of the magnification of ~ 9 000 randomly picked 
source positions at z — 1.5 in the GOODS-fields to a distribution simulated with the SNOC 
package (Goobar et al. 2002) used previously to estimate SN lensing uncertainty in e.g. Knop 
et al. (2003). The distributions, computed assuming NFW halos, are presented in Figure 8, 
where we see that the difference between the distributions for CDFS and HDFN is very small 
and that the mean value of the distributions are again slightly less than unity. The agreement 
between the distributions for the GOODS-fields and the simulated distribution is fairly good 
although the simulated distribution, which does not include any large scale structure effects, 
peaks at a slightly higher value than the one for the GOODS-fields. We conclude that the 
distribution of magnifications obtained using the SNOC package is comparable to a real 
distribution and hence realistic. 



5.3. The Magnification of SN 1997ff 

The issue of magnification of the farthest known SN, SN 1997ff, has been addressed by 
several authors in the past (Lewis & Ibata 2001; Riess et al. 2001; Mortsell et al. 2001; Benftez 
et al. 2002). In the analysis of Lewis & Ibata (2001) two galaxies very close to the line-of- 
sight were considered, both residing at z — 0.56 and visible in Figure 2. Lewis and Ibata 
assumed both galaxies to have the same velocity dispersion and calculated the magnification 
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for three different values. For velocity dispersions 100 kms -1 , 200 kms -1 , and 300 kms -1 
they found a magnification of —0.084 mag, —0.38 mag, and —1.16 mag, respectively. The 
velocity dispersions of these galaxies are closer to 100 kms -1 than 200 kms -1 and the most 
likely magnification estimated by Lewis and Ibata is consequently —0.084 mag. Lensing by 
one of these two galaxies was also studied in Riess et al. (2001) concluding a low probability 
for any significant lensing. More foreground galaxies were considered in the estimation of the 
magnification of SN 1997ff by Mortsell et al. (2001), who considered galaxies within 10" from 
the SN. However, due to uncertainties in the Faber- Jackson relation normalization they found 
the range of magnifications to be too large to allow any quantitative statement. Benftez et al. 
(2002) included 6 galaxies within 15" in their analysis and reported a magnification of ~ —0.3 
mag, although overestimating the velocity dispersions. In their re-analysis the magnification 
is significantly reduced (Benftez, N., private communication). However, the above-mentioned 
estimates all use the filled-beam approximation for calculating distances, i.e. the lensing 
galaxies are put on top of a universe with all matter distributed homogeneously. Since in 
this approach, every line-of-sight is overdense compared to a homogeneous universe, \x > 1 
(if the image is primary) and magnifications are overestimated. 

Including galaxies within 60" and using SIS and NFW halo models, we find a magnifi- 
cation of — 0.13^o!o2 niag and -0.18l£c£ mag (68% confidence levels) relative to an homoge- 
neous universe, respectively. 

5.4. Cosmology Fits 

We have computed corrections for 14 of the 157 "gold" SNe in Riess et al. (2004). In 
this section the implications on cosmology fits due to these corrections are investigated. The 
magnitude of SN 1997ff published in Riess et al. (2004) has been corrected for gravitational 
lensing (Riess, A., private communication) and to facilitate comparisons of fits of cosmolog- 
ical parameters with and without corrections for gravitational lensing we have subtracted 
0.34 mag (Benftez et al. 2002) from the magnitude of SN 199711. 

Figure 9 show the results of fits of and to the data in the top panel. In the 
bottom panel the results of fits of Qm and a constant dark energy equation of state wo to 
the data, assuming a flat universe, are presented. Contours are given at 68%, 95%, and 99% 
confidence levels. 

The confidence level contours of the corrected fit in the top panel move ~ 4% along the 
major axis of the confidence level ellipses toward smaller values relative to the uncorrected 
fit. The effects of the corrections are small also on the f2 M and w fit. In both cases, fits to 
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the corrected data are in better agreement with the concordance model than the uncorrected 
data. Noting that the difference between fits to the corrected and uncorrected data is small, 
one should bear in mind that corrections have been applied only to 14 out of the 157 SNe 
and the effects could thus potentially be larger. No uncertainties in the corrections have 
been taken into account in the fits. Neglecting possible selection effects, simulations of the 
entire "gold" sample indicate that the added uncertainty due to lensing on the estimates of 
dark energy is cr^ A ~ 0.04. Assuming a flat universe, the lensing uncertainty on the equation 
of state parameter becomes a Wo « 0.09. 

6. SUMMARY AND DISCUSSION 

In Gunnarsson et al. (2005) a technique for correcting observed magnitudes of point 
sources for magnification by gravitational lensing was presented and it was also shown that 
the scatter due to gravitational lensing can be reduced by this technique. The technique 
has been applied to the published sample of Hubble Space Telescope discovered SNe in the 
GOODS-fields (Riess et al. 2004; Strolger et al. 2004). Our study shows clear evidence 
for magnified and de-magnified SNe. We show explicitly that it is possible to correct for 
gravitational lensing of SNe, thereby decreasing some of the induced smearing in the Hubble 
diagram at high redshift. We find that the mean magnification factor for the 33 SNe in the 
GOODS-fields is very close to unity, i.e. we find no signs of magnification selection effects 
on the sample. The scatter is consistent with the results of simulations in e.g. Goobar et al. 
(2002). For the most magnified supernova, SN 1997ff, we find a magnification below —0.25 
mag with 95% confidence level, i.e. smaller than what previous studies of this SN concluded. 

The effect of the corrections due to gravitational lensing on the current cosmology fits 
is small. We have computed the corrections for 14 of the gold SNe in Riess et al. (2004). If 
the corrected magnitudes for these SNe are used and f^M and fitted to the gold sample, 
confidence contours move ~ 4% along the major axis of the confidence level ellipses toward 
smaller values. 

Since magnified SNe appear brighter, we expect a correlation between the residual mag- 
nitudes in the Hubble diagram after subtracting the best-fit cosmology magnitudes and the 
lensing magnification factor. Rank correlation tests for the Type la SNe for which we have 
estimated the magnification factor, /z, show a correlation between \i and the residual magni- 
tudes of the expected order. However, uncertainties in both \i and the observed magnitudes 
are large and they can be described as having zero correlation within the 68% confidence 
level. 
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Although our method takes into account undetected matter, there are a few caveats 
concerning the distribution of matter. If light does not trace matter, the distribution could 
be different from what we have inferred from the observed galaxies. The center of dark 
matter halos could be different from the center of galaxies. Moreover, the existence of dark 
halos or even dark compact objects not associated with galaxies is unlikely, but cannot be 
excluded at the moment. Similarly, we have neglected the possible non-spherical shapes of 
dark matter halos. However, the small differences in the magnification probabilities we found 
between the NFW and SIS halo profiles, along with tests we have done indicating that 0.5" 
offsets in the positions of the lenses have also negligible impact in the results presented here 
indicate that the technique is robust. Thus, we conclude that correcting for gravitational 
lensing is advantageous when using standard candles to determine cosmological distances. 

The authors would like to thank Joakim Edsjo, Daniel Holz, Saul Perlmutter and Jesper 
Sollerman for helpful discussions during the course of the work. AG is a Royal Swedish 
Academy Research Fellow supported by grants from the Swedish Research Council, the 
Knut and Alice Wallenberg Foundation and the Goran Gustafsson Foundation for Research 
in Natural Sciences and Medicine. CG would like to thank the Swedish Research Council 
for financial support. 
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z 



Fig. 1. — Redshift dependence of the smoothness parameter. The solid and dashed line show 
rj(z) vs. redshift computed for the CDFS- and HDFN-field, respectively. 
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Fig. 2. — Lines-of-sight to four of the supernovae most affected by lensing. The SN is indi- 
cated by a cross at (0,0). Masses of foreground galaxies are proportional to the size of the 
plot symbols. Small, medium sized, and large filled circles indicate foreground galaxies with 
masses m 20 o/M Q < 10 11 , 10 11 < m 2 oo/M Q < 10 12 , and iri2oo/M & > 10 12 , respectively. Red- 
shifts of the galaxies are written next to the galaxies and the superscripts indicate whether 
the redshift is spectroscopic (S) or photometric (P). 
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Fig. 3. — Magnification PDFs of four of the supernovae most affected by lensing. Solid and 
dotted lines show PDFs computed assuming SIS and NFW halo profiles, respectively. 
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Type la Supernovae 



Fig. 4. — Magnification of Type la SNe in the GOODS-fields. Solid squares and circles 
indicate magnifications computed assuming SIS or NFW halo profiles, respectively. The 
error-bars show the 95% confidence level. 
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Core-collapse Supernovae 



Fig. 5. — Magnification of core-collapse SNe in the GOODS-fields. Solid squares and circles 
indicate magnifications computed assuming SIS or NFW halo profiles, respectively. The 
error-bars show the 95% confidence level. 
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Fig. 6. — Magnification of the 33 SNe vs. redshift. Filled circles and squares represents 
Type la and core-collapse SNe, respectively. Error-bars show the 95% confidence level. The 
figure also shows simulated distributions of magnifications as a function of redshift at 68% 
(light gray) and 95% (dark gray) confidence level (Gunnarsson et al. 2005). The dashed line 
indicate the most probable simulated magnifications. 
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Fig. 7. — Histogram of the magnification of 33 SN in the GOODS-fields. The shaded his- 
togram show the expected (normalized) distribution of SN magnifications from simulations. 
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Fig. 8. — Comparison between distributions for the GOODS-fields and a simulated distribu- 
tion. Solid and dashed lines indicate the distributions of magnifications obtained for ~ 9, 000 
randomly picked sources, at z — 1.5, in CDFS- and HDFN-field, respectively. The simulated 
distribution, indicated by the dotted line, was obtained using the SNOC package. 
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Fig. 9. — Cosmological parameter fits to "gold" SNe in Riess et al. (2004) with and without 
corrections for gravitational lensing. The figure shows fits of f2 M and Q\ in the top panel 
and of Qm and a constant dark energy equation of state Wo in the bottom panel. Contours 
at 68%, 95%, and 99% confidence level obtained from uncorrected and corrected data are 
indicated by solid and dotted lines, respectively and the best fits are indicated by solid circles 
and squares. Note that corrections have been applied to only 14 out of 157 SNe. 
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Table 1. Results for Type la supernovae. 



SIS NFW 



SN 


2 






M 


68% c.l. 


95% c.l. 






68% c.l. 


95% c.l. 


2002hr 


0.526 


14 





993 


[0.983,1 


000] 


[0.977,1.012] 





995 


[0.984,1 


007] 


[0.977,1.021] 


2003en 


0.54 


22 





993 


[0.985,0 


992] 


[0.981,0.997] 





996 


[0.986,0 


995] 


[0.980,1.000] 


2003be 


0.64 


26 





996 


[0.981,0 


997] 


[0.975,1.003] 


1 


002 


[0.983,1 


002] 


[0.978,1.009] 


20021g 


0.66 


21 





971 


[0.970,0 


977] 


[0.967,0.981] 





969 


[0.969,0 


977] 


[0.966,0.984] 


2003bd 


0.67 


36 





997 


[0.988,1 


004] 


[0.981,1.014] 





995 


[0.985,1 


004] 


[0.977,1.015] 
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0.71 


33 
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[0.959,0 


967] 


[0.955,0.972] 





962 


[0.957,0 


964] 


[0.953,0.970] 


2002kd 


0.735 


47 





998 
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006] 


[0.979,1.015] 


1 


003 


[0.990,1 


014] 


[0.976,1.023] 


2003eu 


0.76 


54 
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[0.955,0 


964] 


[0.952,0.968] 





958 


[0.952,0 


960] 


[0.949,0.965] 


2003eq 


0.839 


55 





968 


[0.957,0 


972] 


[0.950,0.978] 





964 


[0.953,0 


968] 


[0.946,0.978] 


2003eb 


0.899 


55 





979 


[0.970,0 


987] 


[0.963,0.996] 





972 


[0.965,0 


984] 


[0.958,0.994] 


20031v 


0.935 


48 





957 


[0.949,0 


963] 


[0.943,0.998] 





952 


[0.945,0 


958] 


[0.940,0.984] 


2003es 


0.97 
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1 


042 


[1.028,1 


070] 


[1.010,1.094] 


1 


068 


[1.049,1 
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[1.031,1.133] 


2002ga 


0.99 
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1 


004 


[0.988,1 


012] 


[0.976,1.025] 


1 


013 


[0.991,1 


026] 


[0.976,1.043] 


2002ki 


1.141 


52 





952 


[0.939,0 


956] 


[0.932,0.964] 





953 


[0.941,0 


960] 


[0.932,0.968] 


2003az 


1.27 


43 





921 


[0.912,0 


925] 


[0.905,0.933] 
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[0.910,0 


924] 


[0.905,0.933] 
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1.30 


87 





948 


[0.939,0 


955] 


[0.930,0.966] 
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[0.932,0 


953] 


[0.925,0.963] 
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1.305 


65 
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[0.947,0 


972] 


[0.936,0.986] 
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[0.945,0 


978] 


[0.934,0.993] 
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1.34 
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[0.971,0 


999] 
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[0.969,1 


001] 
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1.40 
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930] 


[0.912,0.936] 





922 


[0.915,0 


928] 


[0.910,0.934] 


1997S 


1.755 


93 


1 


127 


[1.059,1 


146] 


[1.025,1.186] 


1 


177 


[1.091,1 


207] 


[1.035,1.259] 



Table 2. Results for core-collapse supernovae. 



SIS NFW 



SN 


z 








68% c.l. 


95% c.l. 






68% c.l. 


95% c.l. 


2003ba 





29 
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994 


[0.990,0.994] 


[0.989,0.997] 
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[0.990,0.995] 


[0.989,0.998] 
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[0.980,0.990] 


2003N 





43 


16 





983 


[0.979,0.984] 
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17 
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[0.978,0.995] 
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[0.978,0.996] 
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[0.973,0.999] 
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66 
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1 


120 
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83 
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[0.946,0.975] 





960 
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89 


48 


1 
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[1.001,1.043] 
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1 
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[0.984,1.082] 
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95 


52 





952 


[0.944,0.957] 


[0.938,0.965] 





948 


[0.941,0.954] 


[0.936,0.961] 



